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Abstract 



The free energy of the Coulomb Gap problem is expanded as a set of Feynman 
diagrams, using the standard diagrammatic methods of perturbation theory. 
The gap in the one-particle density of states due to long-ranged interactions 
corresponds to a renormalization of the two-point vertex function. By collect- 
ing the leading order logarithmic corrections we have derived the standard 
result for the density of states in the critical dimension, d = 1. This method, 
which is shown to be identical to the approach of Thouless, Anderson and 
Palmer |^ to spin glasses, allows us to derive the strong-disorder behaviour of 
the density of states. The use of the renormalization group allows this deriva- 
tion to be extended to all disorders, and the use of an epsilon-expansion allows 
the method to be extended to d = 2 and d = 3. We speculate that the renor- 
malization group equations can also be derived diagrammatically, allowing 
a simple derivation of the crossover behaviour observed in the case of weak 
disorder. 
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I. INTRODUCTION 



The phenomenon of the Coulomb Gap has been known for 20 years and has been observed 
both in disordered semiconductors below the Anderson transition and in the high mag- 
netic field Quantum Hall fluid ||^. In any system of localized electrons, when electrostatic 
interactions are taken into account, a soft gap appears in the density of states at the Fermi 
level. This classical effect is the result of the interplay of disorder, long-ranged interactions, 
and the discreteness of electric charge. The Hamiltonian of the system [|T],0 is 

H[n,] = <P^n, + E ^""^ ~ (1) 

where the rij are occupation numbers. The hamiltonian contains two terms; the first de- 
pendant upon disorder and the second upon Coulomb interactions. The 0j are uncorrelated 
energies representing the effects of disorder, and macroscopic thermodynamic quantities 
must be obtained by averaging over them according to the distribution: 

There are two dimensionless parameters governing the problem. The one describing the 
relative strengths of disorder and interactions is 

where a is the lattice spacing, so that Eq gives a measure of the strength of nearest-neighbour 
interaction. We consider the problem on a lattice, although because the system is dominated 

2 

by long-range behaviour over the order of the correlation length (r ~ ^) this should be irrel- 
evant to the physics of the problem. The standard Coulomb Gap problem in strong disorder 
corresponds to 7 <^ 1, whilst 7^1 corresponds to weak disorder. The dimensionless 
parameter describing thermal effects is 

e = ^ = - (4) 
^ T aT ^ ^ 
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Because the distribution of the 0i is symmetric under 0j — — 0j, the Hamiltonian is sym- 
metric under 0j —>■ —(pi^rii —>■ 1 — rii, so that the chemical potential n = and the total 
number of particles is fixed 

Y.n, = N/2 (5) 

i 

We would like to derive the one-particle density of states (OPDOS), g{e), defined as the 
distribution function of the single site energies 

^^.-^E^ (a) 

At nonzero temperatures the problem was treated by Mogilyanskii and Raikh, and in 
particular it was shown that the distribution of electrons over sites obeys Fermi-Dirac statis- 
tics 

^"^^ = ""^'^'^^ = l + exp../T 

Once the OPDOS is known, we can derive many useful macroscopic quantities such as 
screening law and conductivity. 

All standard discussions of the Coulomb Gap focus on the strong disorder problem, although 
we will also consider the case of weak disorder. In this case we might expect a crystalline 
groundstate and a hard gap. However, as noted by Efros 0, the freezing temperature of a 
pure ionic crystal with no disorder is numerically smaller by about two orders of magnitude 
than sq. Thus there is a very wide range of temperatures between this freezing tempera- 
ture and for which the weakly disordered Coulomb Gap system behaves like a strongly 
correlated classical fluid. 

In section we demonstrate that the Coulomb Gap partition function can be expanded as 
a series of Feynman diagrams. The general philosophy is to expand the partition function 
lnZ[0j] as a series of diagrams and only at the last moment to perform the averaging 
over disorder (j^). This method is identical to the method used by Thouless, Anderson, 
and Palmer (TAP) p| to study the spin glass problem; we discuss this in section |T|. In 



the next section we give a simple argument for the lack of any low-temperature phase 
transition in the critical dimensionality d = 1. This means that the derivation of the leading 
approximation to the one particle density of states for the Coulomb Gap problem given in 
section |V| does indeed sum the most important set of diagrams, and that any divergences 
in other sets of diagrams must cancel. The standard results in d = 2 and d = 3 can be 
obtained by means of an epsilon-expansion in e = — 1, as explained in section Finally, 
we speculate that renormalization group theory arguments can allow a derivation of the 
weak-disorder limit in which the standard results for the density of states break down. 



II. DIAGRAMMATIC PERTURBATION THEORY FOR THE COULOMB GAP 



The density of states is calculated by means of a diagrammatic perturbation theory. 
Polyakov shows how to convert classical Ising-like models intoa continuous field problem, 
which can be expanded into the Feynman diagrams well-known from quantum field theory, 
and we use a slightly-modified version to take account of disorder. 
The partition function is rewritten using the identity 



^exp 
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(8) 
(9) 



leaving the result in a form in which the summation over [n, = 0, 1] can be easily performed, 
using the identities: 

- = TT rnsh 

2T 



[nil i 

E exp J2 In cosh 2^ 



cosh = exp ^ In cosh 



2T 



T dX: 

to yield the formally exact result 



^ exp ^ In cosh ^ = exp ^ In cosh 



2T 



(10) 
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exp(2^ — 



dXj 



exp-^lncosh(^^^tA) (n) 



Our general approach will now be to perturbatively expand lnZ[(j)i] as a series of connected 
Feynman diagrams, similar to the 'locator' perturbation series of Anderson 0, and only at 
the very last stage to perform the averaging over 0j. This corresponds to the approach of 
Thouless, Anderson and Palmer (TAP) to spin glasses. 

The diagram rules for the series expansion of InZ are as follows, where nfd{(f)) is the Fermi- 
Dirac function: 

• Each vertex i with lines coming out of it gives a factor of 

• Each line from i to j gives a factor of 

,2 



-e 



rijT 



(13) 



• Symmetry factor l/IG], where \G\ is a combinatorial factor equal to the order of the 
symmetry group of the diagram. This standard factor ^ is well known from statistical 
and quantum field theories. 

When it comes to averaging over disorder, we note that the disorder-averaged bare n-point 
vertex function is: 



Integrating this by parts, we can see that it equals 



x=0 



SO that for all non-pathological distributions of disorder (such as for example Gaussian (|^)) 
the bare two-point vertex function is 



(where d is the dimensionahty of space) and higher vertex functions disappear hke powers 
of 5^ 

An cx ij)--' (17) 

so that they are neghgible to our calculations to leading logarithmic order. In the case of a 
'top-hat' distribution of local site energies, as used in numerical simulations, A„ disappear 
like exp — ^ for n > 2. Thus to leading order, only two-point vertices survive after averaging. 
Note that the bare two-point vertex function is proportional to the non-interacting density 
of states 5^0; and that the renormalized two-point vertex function is given, after averaging, 
by 

(sech^l,) = Ta'g (18) 
where g is the density of states at the Fermi level as a function of temperature. 



III. THE RELATIONSHIP BETWEEN THE COULOMB GAP AND THE TAP 

TREATMENT OF SPIN GLASSES 

The idea of performing a diagrammatic expansion of In Z and of only performing the averages 
over disorder at the last moment is identical to the TAP theory of spin glasses The 
Coulomb Gap problem is in some senses more simple than the spin glass problem, because 
for A ^ T only two-point vertices survive after averaging. The only reason that the field 
theory is not a trivial 0^ theory is the possibility that more than one vertex could represent 
the same site, so that if for example i = k, (pi and (pk will no longer be independant. We 
represent this case by Feynman diagrams including 'connections' (see figure |1|). As in the 
diagrams for the conductivity of a dirty metal |13[, it is these connections which make the 



theory non-trivial, in that there are one-particle irreducible diagrams ||^ more complex than 
a single bare two-point vertex. 
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In the next section we will show that a certain set of Feynman diagrams (see figure ^ which 
we will call the 'maximally crossed' diagrams, are larger by O(^) than any others, and thus 
dominate the expression for the density of states. 

This behaviour is identical to that of spin glasses, for which a certain set of diagrams is larger 
than any others by 0{z) where z, the coordination number, is assumed large. However in 
spin glasses, the sum of a set of loop diagrams incorporating four-point vertices, diverges 
at the spin glass temperature, signalling the onset of replica symmetry breaking (RSB). 
However the TAP approach allows expansions around frozen-in spins, and the method is 
'fail-safe' in the sense that if we are expanding about a state which is not a global minimum 
of energy, we will get negative or zero eigenvectors of the response matrix, and a negative 
or infinite susceptibility. The low-temperature state is marginally stable [|lT| (meaning that 
the eigenvalues of the stability matrix are distributed all the way down to zero) but the 
TAP method holds up under these circumstances. The marginal stability of the system is 
precisely due to the exponentially large number of metastable states. 

In section |V] we note that there exists a simple physical argument why no such phase 
transition can occur at low temperature in the Coulomb Gap problem in d = 1. Thus 
any sets of diagrams whose sums do diverge at low temperatures must cancel with other 
such diagrams. In d = 2,d = 3 the question of whether there exists a glass transition is 
inconclusive, although this does not affect our arguments. 

Anderson demonstrates a derivation of the low-temperature distribution of effective 
local fields in the spin glass problem which is exactly analogous to the standard derivation 
of the density of states in the Coulomb Gap p|. However because of the onset of RSB at low 



temperatures, a rather intricate method due to Bray and Moore |jTT[ is necessary in order 
to perform a perturbation theory expansion around the non-trivial ground-state of the spin 
glass, although their method is a locator perturbation theory, essentially the same as our 
calculation in this paper. Bray and Moore define a Greens function, which when averaged 
over disorder, gives exactly 



The Coulomb Gap is effectively a variant of the random- field ferromagnet with long-ranged 



interactions. Bray [|T2| confirms that whilst the random-field ferromagnet displays low- 
temperature anomalies in the specific heat (again due to marginal stability), it does not 
have RSB. 



IV. ABSENCE OF A LOW-TEMPERATURE PHASE TRANSITION IN D = 1 

As mentioned in the previous section, we are performing a perturbation expansion around a 
trivial groundstate. Thus it is essential that we sum the diagrams at a temperature higher 
than that of any phase transition, and we will show that in d = 1 there is no such transition. 
Consider a 1-d Coulomb Gap system with small but finite disorder, 7 ^ 1, at T = 0. The 
system will form into domains of perfect ionic crystals separated by domain walls. Suppose 
that the density of domain walls is A^, so that the mean number of sites in a single domain is 
A^~^: the limit 7^1 corresponds to ^ 1. Then it is possible to make a simple argument 
similar to that of Imry and Ma , showing that the free energy has a minimum for a finite 
domain wall density, of order 

^2 

AT, ~ (19) 

^0 

and so any disorder, no matter how small, destroys long range order in d = 1. Thus in this 
critical dimensionality there is no low temperature phase transition, and we are safe in the 
knowledge that our perturbation theory is around a stable minimum. 

Vojta [0 shows that within the uncontrolled approximation of the spherical model, the 
Coulomb Gap system has a low-temperature transition to a phase with long-range order in 
(i > 4, whereas in c? < 4 there is no such transition. Monte Carlo simulations are inconclusive 



on the question of whether there is a phase transition in d = 2,d = 3 [T^JT^. Studies of 



a related system, the Coulomb Glass show a phase transition in ci = 3 [13] but not in 



d = 2 [113 . All numerical simulations show clearly that the Coulomb Gap is well-developed 
at temperatures well above the suspected glass transition. This gives us a justification for 
disregarding the glass transition and performing an expansion around the trivial groundstate. 
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V. DERIVING THE DENSITY OF STATES WITHIN PERTURBATION THEORY 



We discuss first the critical dimensionality d = 1, because the divergences encountered here 
are logarithmic and hence tractable to the renormalization group. We mention later how 
the results are modified for the more severe divergences in d = 2,d = 3. We aim to derive 
the results of Raikh and Efros for the Coulomb gap in strong disorder 

9 = 9oo= . , 2 ^° . ,rj.. (20) 

and to show how these are modified in weak disorder. In d=2, where the strong disorder 
density of states is given by 

2T 

= ^ (21) 

Pikus and Efros observed a crossover behaviour for the weak disorder problem, with a 
density of states given by 

A 1 

9 ~ -9oo + ^ exp{-aeo/T) (22) 

We note that heuristically, we expect the low-disorder density of states to behave like the 
Boltzmann function exp — ^ where the typical energy of excitations is aeo. 
We have already discussed how the diagrammatic perturbation expansion of g can be ob- 
tained. The 'maximally crossed' diagrams (see figure for maximally crossed diagrams 
to order Qq ) are larger by O(^) than any diagrams containing unpaired vertices. Before 
averaging over disorder these diagrams give 

(j)' E ■ £ ■ £ ■ f \^ 

sech^(^) = sech^(^) — sech^(^)sech^(^' 



- =ech*(^)=ech*(^) (--J)' - =ech«(^)=ech'(|,, [-^)' + • • • (23) 

We first perform the spatial integration, then by inspection sum the series. As the final step 
of the calculation, we can now perform the average over disorder (^ to give, using (|TB|), and 
to logarithmic accuracy, 
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^-^?o-^ln(l + ||) + --- 
-^?o-e\^ln| (24) 

As we would expect, there is the same logarithmic behaviour for a ^ potential in d spatial 
dimensions. For the Coulomb Gap in d = 1 the divergences are mild enough that we can 
collect together the leading order logarithmic terms depicted in figure ^ to give an expansion 
for g in terms of Qq 

9^ 9o-e go\n— + e gain — (25) 

so that the strong disorder case can be derived as a resummation of all the leading order 
logarithmic terms. A more satisfactory mechanism for collecting together all terms of the 
same magnitude is to use the renormalization group equations, and instead of considering a 
perturbation expansion for g, to derive an expansion the beta-function. This is discussed in 
the next section. 



VI. RENORMALIZATION GROUP THEORY 

We begin the renormalization group treatment by noting that we wish to express the renor- 
malization group equations in terms of dimensionless quantities, such as the partition func- 
tion Z (!]), and 

r = e''a''-^g (26) 

which is the probability that any site is within the gap. We will express T in terms of ^, 
where 

(27) 



where the correlation length associated with a temperature T is p2 

o2 



rT = j (28) 
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The beta-function is 

P{r) = ^ = Po + f3iT + (32T' + --- (29) 

For the case of strong disorder, solving the following equation for the beta-function: 

P{T) = -r (30) 

together with the boundary condition F — > Fq as ^ gives 

To 



F 



1 + Folne 



- l + e4°olnf ^^^^ 
which is exactly the result of Raikh and Efros for the density of states in d = 1 with strong 
disorder ([20|). In other dimensionalities, an epsilon-expansion gives the beta-function to be 



/?(F) = -(rf-l)-F (32) 
which can again be solved with the boundary condition F ^ Fq as ^ ^ to give 



F ocf 



-(d-i) 



9t{0) oc ^ (33) 

which is the standard result of Efros and Shklovskii 0]. Thus we reach our main conclusion, 
that the TAP method for solving the spin-glass problem is an identical approximation to 
the Efros-Shklovskii method for obtaining the density of states of the Coulomb Gap. 
We note that the d — 1 term can be derived using hyperscaling. The system has a critical 
point which has been moved to T = 0, and the correlation length (p8|) implies a critical 
index = 1, so by hyperscaling, the heat capacity must be proportional to T*^"^ and the 
density of states to T'^"^. 

In the limit of weak disorder 7^1, the standard results for the density of states break down 
2T| , and in this section we explain how to use the renormalization group equation to derive 



the behaviour in this regime. The renormalization which we are performing corresponds to 
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a change in the disorder, 7 of the system. As noted in section 0, we expect a Boltzmann-hke 
behaviour for the density of states at low disorder. 

r = r* exp ( — ^ I 



T J 

inr/r, = 

/3(r) = inr/r, (34) 

where F^, is an order unity constant. This beta-function is a result which we expect to be 
universal at low disorder. By solving equation (^) for the modified beta-function which 
crosses over between the weak-disorder and strong-disorder limits, 

/3(r) ^ -(rf- 1) -r (r < 1) (35) 
/3(r)^inr/r, (r~r,) (36) 

with again the same boundary conditions F — > Fq as ^ 0, we obtain a density of states 
with a crossover behaviour which md = 2 coincides with that (^) observed in the numerical 
simulations by Efros and Pikus [ PT[] . 

All the information in equations can be summarized in figure ^.Any glassy phase 

transition or crystallization such as those discussed in section |IV| would constitute a breaking 
of universality in the region InF ~ 0. However even for infinitesimal disorder, the site 
energies will have a spread of order Eq and F will be of order unity so that this region is of 
little physical significance. 

As noted earlier, there is a crossover between the strong-disorder (Efros-Shklovskii) and 
weak-disorder (Boltzmann) behaviour. In d = 1 this occurs at a value of 

r = F, = a + 61nF, (37) 

with a, h order unity constants, whilst md = 2 the modified renormalization group equations 



give exactly the DOS observed in [21 



The corrections to hyperscaling in the RHS of ( |2^ ) can be obtained as a systematic pertur- 
bation series by using the identity 
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(i(lnZ) 
d\nT 




(38) 



which follows from ([TT|). Thus we can develop a perturbation theory for (3, the first non-zero 



term of which (see figure |^) corresponds to the standard Coulomb gap result (|30| , |35|) . We 
believe that by evaluating further terms in the series that we will obtain a series which 
interpolates smoothly between ( [55| ) and (PB]). 



The essence of our paper is in linking two separate pieces of physics, the spin glass problem, 
and that of the Coulomb Gap. We have demonstrated that the standard results for the den- 
sity of states in the Coulomb Gap can be derived by a method identical to the TAP method 
in spin glass theory.The Coulomb Gap density of states plays the part of the renormalized 
two-point vertex function in an effective field theory. In = 1 we can perform the renormal- 
ization directly by summing a dominant set of diagrams to obtain the strong disorder limit 
of the Coulomb Gap. The results in d = 1 suggest the use of the Renormalization Group 
to derive this result in any dimensionality, and to extend the result into the case of weak 
disorder. It is shown how a diagrammatic expansion for the /9-function can be used to derive 
the renormalization group equations. Thus we have demonstrated that our formalism allows 
an understanding of the cross-over between the strongly disordered Coulomb Gap problem 
and systems with weak disorder. 
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FIGURES 

FIG. 1. An example of a Feynman diagram containing a connection; this represents the possi- 

biUty i = k 



FIG. 2. Diagrams giving to the standard Coulomb Gap density of states; these 'maximally 
crossed' diagrams are larger by 0{j-) than all others. 



FIG. 3. The Dyson equation giving all the leading logarithm terms in the density of states 



FIG. 4. The scaling beta- function ci In F/ci In ^ plotted as a function of InF 



FIG. 5. The first non-zero term in the perturbative expansion of the beta-function. The wavy 
line represents a ^d/d^ term 
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